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Abstract 

Three forms of the high-frequency asymptotic Green's function for Lilley’s equation are reviewed and 
compared to the exact solution over a wide range of Strouhal numbers. The asymmetric approximation, 
which applies to sources away form the jet axis, and the quasi-symmetric approximation, which is arrived 
at by making a near-axis source assumption, are both obtained for parallel round jets from a formal 
Fourier-transform solution. The ray-theory solution, which is the only high-frequency approximation that 
can be applied to more general mean flows, follows from a WKB ansatz and is shown to be closely related 
to the asymmetric approximation. The comparisons show that the best overall prediction of the exact 
Green’s function is given by the asymmetric approximation which remains accurate down to a Strouhal 
number of 1 /2. The close relationship between the asymmetric and ray-theory approximations suggests 
that the high-frequency asymptotic Green’s function for more general mean flows would be similarly 
successful. 

1 Introduction 

The prediction of jet noise has been an area of continuous interest over the last half century and has become increas- 
ingly important in recent years due to stricter noise regulations placed on the commercial aircraft industry. The need 
for quieter jet engines has led to several noise abatement techniques such as the placement of tabs and chevrons at the 
nozzle exit in order to alter aerodynamic sound generation through enhanced mixing. Accurate and robust prediction 
tools are instrumental in the design of more efficient noise- suppression devices. 

Noise generation in jets is generally regarded as a by-product of the unsteady features of the flow and, in many sit- 
uations of practical interest, the dominant feature is turbulent mixing. Jet mixing noise emanates from both fine-scale 
turbulence and the unsteady motions of large-scale coherent structures. The contribution from the latter component is 
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usually most important at shallow angles off the downstream axis (especially in supersonic jets) and is often success- 
fully predicted using either an instability-wave based approach or a large eddy simulation. It is this noise source that is 
most directly impacted by mixing enhancement devices since they tend to breakup the large structures at the expense 
of creating more fine-scale turbulence. 

The present investigation is concerned with the mixing noise due to fine-scale turbulence which dominates the 
spectra away from the downstream jet axis. This component is most often analyzed by employing the acoustic analogy 
and assuming that the noise generating eddies are compact and behave as convected acoustic sources. The sound 
field is then governed by Lilley's equation which describes the acoustic propagation on a specified mean flow due 
to multipole-type sources. The source distribution is, in general, modelled using appropriate space-time correlation 
functions. 

Solutions to Lilley's equation are typically constructed by introducing a Green's function. This allows the mean- 
flow refraction effects to be determined independent of the source distribution and limits the empiricism inherent 
in the acoustic analogy to the source modelling problem. Despite the linear nature of Lilley's equation, accurate 
numerical determination of the associated Green’s function for an arbitrary mean flow is still a major undertaking 
[1] and consequently much attention has been focused on the simplifications to be gained by use of high-frequency 
asymptotics [2, 3,4]. 

The high-frequency limit arises when the acoustic wavelength of the aerodynamic noise is much shorter than the 
characteristic length scale of the mean flow. For simple round jets, this is usually the case within the first several jet 
diameters downstream of the nozzle exit where the mixing layer is thin and the turbulence intensity is at its peak. 
The high-frequency noise can be further increased by the presence of mixing enhancement devices. As a result, 
high-frequency asymptotic approximations to the Green's function for Lilley's equation form the backbone of many 
jet-noise prediction schemes, eg. the so-called MGB computer code [5] and its derivatives. 

For arbitrary mean flows, the high-frequency solution to Lilley's equation is described in terms of the ray-theory 
of acoustics [4]. However, it is often reasonable to assume, for high Reynolds-number jets of practical interest, 
that the mean flow is both locally parallel and axisymmetric - even jets issuing from tab and chevron nozzles are 
known to become axisymmetric within 4 to 5 diameters downstream of the nozzle exit. When the locally parallel and 
axisymmetric assumption is made, three different closed-form expressions for the high-frequency Green's function for 
Lilley's equation are found in the literature. Goldstein [6] developed an approximation for the Green's function by 
restricting attention to sources located several acoustic wavelengths off the jet axis and determining the high-frequency 
asymptotics of the formal Fourier-transform solution available for parallel round jets. Balsa (2, 5] also obtained an 
expression for the high-frequency Green's function from the formal Fourier-transform solution but did so by assuming 
that the source lies near the jet centerline. Finally, Goldstein [3] presented a closed-form ray-theory solution for 
parallel round jets. 

The goal of limiting the empiricism in jet-noise prediction schemes to the source modelling problem will be 
achieved only when the high-frequency asymptotics provide an accurate approximation of the exact Lilley's equation 
Green's function. Adequate agreement between the exact source directivity and a high-frequency approximation for 
Strouhal numbers as small as one has been demonstrated in some limited circumstances by Tester and Morfey [7] using 
a ray-theory solution for polar angles outside the zone of silence of a round jet and by Scott [8] for a two-dimensional 
isothermal flow with a piecewise constant mean shear. The primary objective of the present paper is to determine the 
relative success of the above three high-frequency approximations for parallel round jets by comparing them to the 
exact order-one frequency solution over a wide range Strouhal numbers and far-field observation angles. In doing so, 
the relationships between the different solution forms will be revealed and the potential success of the high-frequency 
approximation for more general mean flows will be discussed. 

The general problem defining the Lilley’s equation Green's function for a uni-directional transversely sheared 
mean flow is presented in section 2 where the formal Fourier-transform solution available for parallel round jets is 
given. The high-frequency asymptotic behavior of that formal solution is considered in section 3 where the approxi- 
mations corresponding to the analyses of Goldstein [6] and Balsa [2, 5] are summarized. Appendices A and B provide 
the details of those analyses with the former correcting an error in the derivation of Goldstein [6]. A comparison of 
the two high-frequency approximations with the exact order-one frequency Green’s function is given in section 4. It 
is shown there that the best over all prediction of the exact result is provided by the corrected expression of Goldstein 
which is referred to here as the asymmetric high-frequency approximation and which remains accurate down to a 
Strouhal number of 1/2. The Balsa expression is referred to as the quasi -symmetric high-frequency approximation 
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and is shown to be at its best when applied to the ring-source directivity in which case it becomes coincident with the 
asymmetric approximation at sufficiently large polar angles off the downstream jet axis. 

The ray-theory solution, which is reviewed for a uni -directional transversely sheared mean flow in Appendix C, is 
shown to be closely related to the asymmetric high-frequency approximation in section 5 where the generalization to 
complex rays in the zone of silence is also considered. In section 6, local modifications to the ray-theory solution near 
the caustic and branch point are constructed and a composite solution is presented. The composite ray-theory solution 
is shown to provide a good approximation to the asymmetric high-frequency solution over the entire range of Strouhal 
numbers considered and from this result it is inferred that the high-frequency approximation for more general mean 
flows (i.e. mean flows that are non-axisymmetric and/or non-parallel) would be reasonably accurate down to Strouhal 
numbers as small as 1/2. 

2 Formulation 

Interest here is in the acoustic propagation on a parallel, doubly infinite jet for which the mean flow is given by 

u = iu(y,z), P — p(y ?2), c = c(y. z), p = constant. (1) 

where u y p , c and p are the velocity, density, sound speed and pressure, respectively and (1) is an exact solution to 
the inviscid, non-heat-conducting equations of motion. The Cartesian coordinates x = {x.y.z} are chosen such that 
x is aligned with the direction of the mean flow and the unit vector i is in that direction. The mean-flow profiles are 
required to approach constant ambient values, 

U->0, C-*Coo, 


as yjy 1 + z 2 — * oo. 

Assuming a calorically perfect ideal gas, the linearized equation governing the acoustic propagation on (1) is [9] 

“ = § i{w^ v -^) n+2elvi - v i n = r ’ ® 

where 77 denotes the acoustic pressure fluctuation normalized by pc 2 , 

D_ = d_ _d_ 

T>t ~ dt^ U dx 

is the convective derivative relative to the mean flow and t denotes the time. The term r represents the acoustic source 
distribution and is given by 

r = h v - f - 2 ™&' (3) 

when this quantity is produced by a fluctuating force per unit volume. In the absence of temperature fluctuations, 
Lilley’s equation [10] is obtained by replacing / with the quadrupole source distribution f = V*(v <g> v) where v is 
the velocity fluctuation relative to the mean flow and ® denotes the tensor product. 

Since (2) is linear, the solution for an arbitrary source distribution can be obtained through superposition of 
solutions to 

£j {Gu{x\x s )e~ lu>t ] = [c^,5(x - x s )e -,u, ‘] , (4) 

where u) is the frequency, x 3 is the source position, S is the Dirac delta function and G u denotes a reduced Green’s 
function. It is common practice to include the convective derivative D /D t in the inhomogeneous term of the reduced 
Green’s function equation because doing so simplifies the subsequent computation of the acoustic field when attention 
is restricted to the first term in (3) - the so-called self-noise term. It should be noted however that use of (4) does not 
limit the form of r since the reduced Green’s function G^ corresponding to the right-hand side, 

clo5(x - x s )e~' ut , 
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is related to G w by 


(5) 


(a? ]*c$) — -j - u s ^ ^u> I *^$) > 

where the subscript 5 denotes evaluation at the source position. 

When the mean flow depends only on the radial coordinate in the y-z plane r = yjy 2 + z 2 , a formal solution for 
can be obtained by reducing (4) to a system of linear ordinary differential equations. Following Goldstein [6], the 
reduced Green's function is written as 


1 +OC p + oc 

G u (x\x s ) = j^ £ e ^v-V') j Gnirl^^k^c-^-^dkx, 


( 6 ) 


where <p = arctan (z/y) is the azimuthal angle in the y-z plane. The Fourier coefficients G n are determined by 


<P 2 d 

r dr 




6(r - r s ) 
rn 2 


(7) 


where 0 = (l + KM)/a,ko = ^/coc, k = fci/fco, a = c/cqo is the local sound speed normalized by its ambient value 
and M = u/c ^ is the local Mach number based on the ambient speed of sound. 

Equation (7) must be solved subject to the conditions that G n remains bounded at r = 0 and behaves like 
an outgoing wave as r — > oo. The solution satisfying these conditions can be expressed in terms of two linearly 
independent homogeneous solutions to (7), say w\ and w 2 , as follows [11] 


Gni^r^uj.kx) 


W] ( r\n)w 2 (r a |/c) 
r s a 2 s W{r s |«) 


for r > r s . 


( 8 ) 


where 

W(v\k) = wi(r\tx)w f 2 (r\K) — w[(r\n)w 2 (r\K, ) 

is the Wronskian, a prime denotes differentiation with respect to r and the w 3 have been chosen such that 


w i 


► constant x r K 

w 2 — ► constant x r' n ', 


. as r - 
as r — > 0. 


oo, 


(9) 

( 10 ) 


Only the r > r 3 form of G n is given since primary interest is in the behavior of the solution in the far field. 

The homogeneous solutions w 3 must be determined numerically in general. However, when the frequency is 
sufficiently large, the equations become simple enough to be solved analytically. The high-frequency limit has the 
additional benefit of ‘cutting-off’ any spatially growing instability waves which appear as discrete eigensolutions to 

(7). 


3 High-frequency, far-field approximation 

The high-frequency limit describes the situation wherein the wavelength of the acoustic field, 1/ko ~ c 0 0 fu>, is 
much shorter than the characteristic length scale of the mean flow, viz. the jet radius 



HW d r = — 

c(r) Mj 



M{r ) 
a(r) 


dr, 


(ID 


where Mj is a jet Mach number which is taken here to be the ratio of the axial velocity to the local speed of sound at 
the jet centerline. Since the mean flow is parallel in the present analysis, the streamwise wavelength l/k\ scales like 
l/fco and the high-frequency limit can therefore be expressed mathematically as 


r— = ko oc. with — = k = 0(1), 

Cqq k() 

where it has been assumed, for simplicity, that the mean-flow quantities M, a and rj are all order one. 
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In many technological applications, one is only interested in the behavior of the acoustic field at remote distances 
(in terms of the characteristic mean-flow length scale) from the aerodynamic noise sources. The so-called far-field 
behavior is most conveniently expressed by introducing polar coordinates in the x-r plane, 

R = yj(x- x s ) 2 + (r - r 3 ) 2 , 9 - arccos ? 

with the origin at the source point, and considering the limit as R — > oc. 

Using (9), it can be shown that the integrand in (6) has a point of stationary phase at 

k = —cos 9 + 0(R~ y ) : 


as R — ► oo. The integral over k\ can then be approximated using the method of stationary phase [11] with the result 
that 



G n (r|r*;w,fci)e lkl(x Xa) dki 


( 2^-fc 0 sin 2 5 wi(r| - cos 0)w 2 (r s \ - cos 6) iko R cos 2 0 
\ i R ) r s a 2 W(r s \ - cos 9) 


( 12 ) 


as R ^ oo. 

The asymptotic approximation for is completed by determining the high-frequency behavior of W\ and W 2 . 
To facilitate this, new dependent variables vi and V 2 are introduced as follows 

$(r) 

Wj ( r | — cos 9) — ^2 v 3 (r) for j = 1,2, 

where now # = (1 — Mcos9)/a. Substituting into (7), (9) and (10) shows that the v 3 must satisfy 

v" H- (k$Qn + S)t> = 0 

subject to 

v\ — ♦ constant x e lfcorsin0 , as r — ► oo, 

V 2 — ► constant x as r — ► 0, 

where 


rQ n (r) = \/r 2 q 2 - Jn/ko) 2 , q(r) = V ~ cos 2 0, 
x _ $ ( r<P'\ f ( 1 \ 2 
s(r) = ) + (^) ! 


(13) 

(14) 

( 15 ) 

(16) 
(17) 


and the square roots in ( 1 6) are chosen such that they have positive imaginary parts for negative arguments. Introducing 
v\ and V 2 into (12) and the result into (6) leads to 


G UJ (x|x s ) 


as fco, R — * ► oo, where 


G<jJ I 'Es')R'UJ ) 
a s (l—M s cos9) ' 

pikoR 


^ (x|Xs)_ 4nR ’ 


(18) 


(19) 


corresponds to the reduced free-space Green's function when both the source and far-field observation points lie in the 
same azimuthal plane, 




mAip—ikoRsin 2 9 


( 20 ) 
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Aip = ip — (f s and V = viv f 2 — v[v 2 is the Wronskian which, in view of (13), is independent of r. 

Before the high-frequency asymptotic solutions to ( 13) can be constructed, the scaling of the azimuthal wavenum- 
ber n with ko must be considered. The summation in (20) suggests the need for uniformly valid approximations to v\ 
and v 2 for all n. However, the dominant behavior of in the limit as ko — ► oo is determined by a relatively small 
number of azimuthal modes centered about a critical value of n and it suffices to construct asymptotic solutions to ( 1 3) 
based on the k 0 scaling of that critical value. 

When the distance between the source and the jet axis is sufficiently large (i.e. several factors of I/Ajq), the result- 
ing acoustic field is asymmetric and the critical azimuthal wavenumber behaves like the streamwise wavenumber and 
scales with ko . This scaling was considered by Goldstein [6] and will be referred to as the asymmetric, high-frequency 
approximation. As the source moves closer to the jet axis the acoustic field becomes increasingly axisymmetric and 
eventually the critical value of n scales like l/r j. The near-axis source problem was analyzed by Balsa [2, 5] and will 
be referred to here as the quasi -symmetric, high-frequency approximation. 

Solutions to (13) using the asymmetric, high-frequency scaling, 

n = O(ko) as ko — ► oo, 


were constructed by Goldstein [6] however the bounded solution used in that analysis (equation (6.61) of reference 
[6]) is in error. The correct asymptotic solutions for v\ and V 2 are given in Appendix A. Using those results in (20) 
leads to 


n u 


-hoo 

E 


2 yZ-Vnirs) 
ko ^sQni^s) 


i 


Ai[7?„(r 5 )]e 


mAv?+ifco(Cn — -ftsin 2 &) 


(21) 


as ko, R — i ► oo, where 


Vn(r) = - [ffcoCn(r)] * , C n(r) = [ Qn(r)dr. Ql(n) = 0, 

J r$ 

Ai denotes the Airy function and the cube root in the definition of r) n is taken such that rj n ^ 0 for r ^ r$. 
Solutions to (13) using the quasi-symmetric, high-frequency scaling, 

n = 0(1) as ko —> 1 oo, 

are given in Appendix B. Substituting (59) into (20) yields 

+oo 


( 22 ) 


as ko, R —> oo, where 

S(r)= [ q(r)dr = Co (r) ~ Co(0). 

Jo 

H^ i:) and 3 n denote the Hankel and Bessel functions of the first kind, respectively, and the quantity under the square 
root in (22) has been simplified by making the approximation « r s q s which was also used in reference [5] and 
becomes increasingly accurate as r s — ► 0. 

One of the advantages of the quasi-symmetric approximation is that the sum over n in (22) can be evaluated in 
closed form. Using Graf’s addition theorem [12], one can write 

+ OC 

Y, a^(koOin(ko^y nA ^ = Hq 1 * (^v'^ + €?-2€6cosA V >) , 


which, when substituted into (22), leads to 


^ e ifc 0 (£-Hsin 2 &~€s cos A <p) 

as ko, R — ^ > oo, where H[ ) 1) has been replaced by its large argument behavior since £ — » oo in the limit R —> oo. 


(23) 


NAS A/CR— 2003-2 1 2089 


6 



Figure 1 . Mean-flow profiles (25) and (28) with Mj = 0.9 and Tr = 3. 


4 Comparison of exact and asymptotic solutions 

Comparisons of the high-frequency, far-held approximations to the exact solution for the reduced Green's function 
are most easily carried out in terms of G w . It follows from (5) and ( 1 8) that 




i G u (a? | tCj ) 

<jo( 1 — M s cos9) 


(a? |a?g (a? | x $ ) 

iva s ( 1 - M s cos9 ) 2 


(24) 


as ko , R — ► oo, where is given by (2 1 ) or (23) for the asymmetric or quasi-symmetric approximation, respectively. 

The exact solution for is obtained numerically using the adjoint Green's function scheme given by Tam and 
Auriault [13]. The interested reader is referred to that reference for the details. 

The mean-flow Mach number profile used in the comparisons is 


u(r ) 
c(r) 


M(r) 

a(r) 


Mjsech 2 ( 2 r), 


(25) 


and it follows from (11) that the jet radius r j — 1/2. The profile is shown in figure 1 for a jet Mach number M j = 0.9. 


Two different profiles are considered for the sound speed ratio a. The first corresponds to an isothermal jet which, 
in view of the ideal-gas result 

c = v / 7 (26) 

has a constant speed of sound and therefore 

a(r) = 1 , (27) 

where 7 is the ratio of specific heats, 5R is the gas constant and T is the mean static temperature. The second profile is 
obtained from (26) and the Crocco-Busemann law and is given by 


where 


a 2 (r) = 1+11 + 


2 — -Mj) 2 — J_ 
2 J ) MjJTr 


M(r) - 2-±M 2 (r), 


Tr = 1 + 


7“ 


b 


(28) 
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Figure 2. Turning point functions for an isothermal jet at various 6 . 


is the ratio of the stagnation temperature at the jet centerline to the ambient temperature- In the results presented here, 
7 = 1 . 4 , = 287. 06 J /kgK and = 290K. The a profile obtained from (25) and (28) with Mj = 0.9 and Tr — 3 

is shown in figure 1. 

In deriving the asymmetric, high-frequency approximation (21), it was assumed (see Appendix A) that, for each 
order- k 0 value of n, (13) has one simple turning point rs corresponding to a zero of Q 2 n . Figure 2 (a) is a plot of r 2 q 2 
for various polar angles 6 using (25) and (27) with Mj = 0.9. Similar curves are obtained for hot jets with a given by 
(28). Since the zeroes of are determined by solutions to 

rV = (n/ko) 2 . (29) 

figure 2(a) shows that the turning-point assumption made in Appendix A holds for 0 < 9 < 150 provided n/ko ^ 0. 
For 0 > 150, multiple turning points are possible for certain values of n/ko. This behavior can be accounted for by 
making appropriate modifications to the analysis in Appendix A, however these upstream angles are usually of no 
practical interest and will not be considered here. 

The curves in figure 2(a) show that the turning-point assumption of Appendix A is violated when n = 0 and 
0 > 58 in which case (29) has no real solution for r$ (see also figure 2b). Since the azimuthal wavenumber scaling 
n = O(k 0 ) clearly does not apply to the n = 0 term in (20), it would seem that this term must be evaluated using 
the n = 0(1) scaling of the quasi- symmetric approximation. However, if the convention that rs = 0 when Q ^ has 
no zeroes is adopted (as is done for r a in Appendix B), it can be shown that (54) and (59) agree at n = 0 by simply 
replacing the functions Ai, Hq 0 and Jo appearing in these expressions by their large argument behaviors. The same 
proof also applies when Ql has a zero provided the simplifying assumption r$ — r s /cq 2/3 of Appendix B is made. 
It therefore follows that the results of Appendix A remain valid even when n = 0. 

The quasi -symmetric analysis of Appendix B is based on the assumption that (56) has at most one simple turning 
point tv corresponding to a zero of q 2 . Figure 2(b) is a plot of q 2 for various angles 0 and the same isothermal mean- 
flow profiles used in part (a). Again, similar curves were found for heated jets. It is clear from the figure that q 2 has 
one simple zero for 0 < 0 < 58 and none for 6 > 58. At 6 « 58, q 2 has a higher-order zero at r = 0 which was not 
accounted for in the analysis of Appendix B. The effect of this omission is highly localized and, apart from a small 
neighborhood of 9 « 58, the turning point assumption of Appendix B is satisfied for the entire range of polar angles 
considered here. 

When comparing the exact and asymptotic solutions for G w , it is convenient to work in terms of the ratio G u /Gw 
since this quantity becomes independent of the radial parameter R in the high-frequency, far-field limit. Figures 3-6 
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Figure 3. jG^/^l x 10 4 for an isothermal jet with r s = 0.75 and St = 2. Solid lines, asymmetric approximation; dashed lines, 
quasi-symmetric approximation; symbols, exact solution, (a) A = 0; ( b ) A<p = 30, dot-dashed line, composite ray solution; (c) 

A <p = 90; (d) A<p = 120. 


show \Qtj/Guj \ as a function of the polar angle 0 for various values of the Strouhal number, 

st = u 2rj = k ° Tj 

2n Cooa(0)Mj n a(0)Mj 

Parts (( a-d) of each figure correspond to azimuthal-angle parameters Atp of 0, 30, 90 and 120. The asymmetric and 
quasi-symmetric approximations are indicated by the solid and dashed lines, respectively, and the symbols correspond 
to the exact solution. The asymmetric approximation given by (21) was evaluated by summing the convergent series 
in the azimuthal wavenumber n from —N to +JV where N varied from 11 at St = 2 down to 4 at St = 0.25. The 
results in figures 3-6 were computed for a point source at r 3 = 0.75 using a mean flow given by (25) and (27) with 
Mj - 0.9. 

Figure 3 shows excellent agreement between the asymmetric approximation and the exact result for all A<p. The 
level of agreement tends to diminish as St decreases but remains fairly good down to St = 0.5, cf. figure 5. At 
the smallest Strouhal number (cf. figure 6), the asymmetric approximation, although no longer in good quantitative 
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Figure 4. As figure 3 but for St = 1. 


agreement, does still correctly predict the trends of the exact solution - a peak near 0 = 60 with a steep drop off for 
0 < 60 and a more gradual decline for 6 > 60. 

Figures 3-6 show that, at 0 = 90, both the asymmetric and quasi- symmetric approximations remain in near perfect 
agreement with the exact solution for all St . At 0 = 90, the mean-flow refraction effects for an isothermal jet vanish as 
can be readily verified from the expression for <P. It is interesting to note that this exceptional case is captured equally 
well by both high-frequency approximations despite their differences in assumed azimuthal wavenumber scaling. 

Considering the quasi -symmetric approximation, the most striking feature revealed by figures 3-6 is the failure 
to predict the oscillatory behavior of the exact solution in the range 30 < 0 < 60. This is due in large part to the 
assumption introduced in reference [2] (and used in Appendix B) that the no-tuming-point form of G w adequately 
approximates the Green's funcdon when r s > r a > 0. It can be seen from figure 2(b) that (56) has a turning point 
that lies between the source position r s = 0.75 and the jet axis when 30 < 9 < 60. By neglecting the turning point, 
the quasi-symmetric approximation does not allow for ray-interference effects which are primarily responsible for the 
oscillations in the exact solution. The issue of interference is discussed in more detail in section 6 where the ray -theory 
approximation of the Green’s function is considered. 

The figures also show that, in the range 6 < 30, the quasi-symmetric approximation tends to over predict |G W jQ w \ 
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Figure 5. As figure 3 but for St = 0.5. 


when A <p < 90 and under predicts this quantity when A > 90. For 9 > 60, the agreement between the quasi- 
symmetric approximation and the exact solution actually improves as St decreases. This surprising behavior, which 
was also noted by Balsa [2], is probably due to the increase in acoustic length scale with decreasing Strouhal number 
which, for a fixed source position, makes the solution appear more axisymmetric. It should also be noted that the 
ray-interference effects near 9 = 60 diminish as St decreases and this too improves the agreement with the exact 
solution. 

When modelling the acoustic-source distribution r for use in a noise prediction scheme, it is usual to assume 
that the sources are sufficiently compact so that only the absolute value of appears in the resulting formulae. 
Nevertheless, it may be of some interest to see how well the high-frequency, far-field approximations predict the phase 
of Gcj/Gv This quantity is plotted in figure 7 at St = 1 for various A <p and the same mean flow used for figures 3- 
6. The curves show that overall the asymmetric approximation is better than the quasi-symmetric approximation at 
predicting the exact result. The quasi-symmetric approximation is at its best when A <p < 90 and 9 > 45. 

In the case of round jets, it is often also assumed that the strength and orientation of the sound sources that make up 
r are independent of the azimuthal angle. The sound field emitted by such a source distribution is then axisymmetric 
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Figure 6. As figure 3 but for St = 0.25. 


and, at fixed axial and radial positions, can be characterized by a ring- source directivity factor which is defined here as 


V = 


/ +7T 
-7T 


Gu 

6 h 
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u> 2 a 2 (l — M s cosOy 
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for the asymmetric approximation, and from (23) that 
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2 d^ 5 


~ 27re“ 2/£oIm ^ a Io(2/coIm^ s ) « 27re- 2fc ° Im( ^“H 


(31) 


for the quasi -symmetric approximation, as /co, i? — > oo, where Io denotes the modified Bessel function of the first kind 
and zero-th order and a subscript a indicates evaluation at the turning point r = r a . The approximation introduced 
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Figure 7. arg /GJ) for an isothermal jet with r s = 0.75 and St = 1. Solid lines, asymmetric approximation; dashed lines, 
quasi-symmetric approximation; symbols, exact solution, (a) = 0; ( b ) — 30; (c) Ap = 90; ( d) Ap = 120. 


on the right-hand side of (31) corresponds to a leading-order composite expansion for the ring-source directivity and 
puts the quasi-symmetric approximation in agreement with the relation used in the so-called MGB computer code of 
reference [5]. 

The ring-source directivity factor is plotted as a function of polar angle 0 in figures 8-1 1 for Strouhal numbers 
St = 2, 1, 0.5 and 0.25. Parts (a) and ( b ) of each figure correspond to the mean flow given by (25) and (27) with 
Mj = 0.9 but differing source positions, r s = 0.5 and 0.75, respectively. Similarly, parts ( c ) and ( d ) show results at 
r 3 = 0.5 and 0.75, respectively, for the heated jet given by (25) and (28) with Mj — 0.9 and Tr ~ 3. 

The figures show that the agreement between the asymmetric approximation and the exact result is excellent at 
St = 2 and remains good down to St = 0.5. Even at St = 0.25, the approximation does a fair job at predicting the 
exact result for 0 > 60 which corresponds to the range where the jet noise is dominated by fine-scale turbulence in 
many technological applications. At the two largest Strouhal numbers, the level of agreement shows no sensitivity to 
changes in the source radius or mean-flow temperature profile. A sensitivity to the source radius becomes apparent at 
the two lowest St but only in the range 6 < 60. This may indicate an increased importance of the n = 0 term in (30) 
as the source moves toward the jet center line. 
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Figure 8. Ring-source directivity £> x 10 6 for isothermal (a, b) and heated (c, d) jets at St = 2. Solid lines, asymmetric approximation; 
dashed lines, quasi-symmetric approximation; symbols, exact solution, (a, c), r s = 0.5; ( b , d), r s = 0.75. 


The agreement between the quasi-symmetric approximation and the exact solution is best when 6 > 60 regardless 
of St. In this range of 6 , the quasi-symmetric and asymmetric ring-source approximations are nearly coincident - a 
result that is not simply fortuitous. By approximating the sum in (30) with a Riemann integral [11], replacing Ai by 
its large (negative) argument behavior and determining the large- ko behavior of the resulting integral, it can be shown 
that the right-hand side of (30) becomes 27 t as k 0 — s ► oc when r $ - rs » k ^ 2/3 which is in exact agreement with the 
quasi-symmetric approximation (31). 

For 9 < 60, the quasi-symmetric approximation is never very good, primarily because of a consistent over 
prediction of the 6 at which V is a maximum. This may indicate the increased importance of asymmetric effects when 
the Green’s function problem possesses a turning point. A conjecture that is supported by the fact that the agreement 
in the downstream range 6 < 90 appears to be worse for the heated jet which, it turns out, has a broader range of 
6 where a turning point arises in the Green's function problem. The improvement in overall agreement between the 
quasi-symmetric approximation and the exact result observed with decreasing St in the point-source comparison is not 
apparent in the ring-source comparison which suggests that the point-source observation may be strongly dependent 
on the azimuthal angle parameter A ip. 
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Figure 9. As figure 8 but for St = 1. 


5 Connection to Ray theory 

The comparisons of the previous section show that, for parallel round jets, the asymmetric, high-frequency ap- 
proximation is in good agreement with the exact Lilley's equation Green’s function over a wide range of Strouhal 
numbers. One would like to demonstrate a similar level of agreement between the exact and high-frequency asymp- 
totic solutions for more general mean flows. Unfortunately, when the mean flow is neither axisymmetric nor parallel, 
the problem governing the acoustic propagation can no longer be reduced to a system of linear ordinary differential 
equations by Fourier analysis and an asymptotic analysis of the type described in section 3 is no longer possible. Nev- 
ertheless, analytic progress can still be made for more general mean flows by considering the high-frequency limit. 
The corresponding asymptotic solutions are then described in terms of the ray theory of acoustics [3, 4]. 

The ray-theory solution of (4) for the uni-directional transversely sheared mean flow (1) is summarized in Ap- 
pendix C. It is implicitly assumed, when constructing this solution, that the length scale of the acoustic field is order 
1 /ko in all directions. Consequently, the ray-theory approximation (77) is expected to be closely related to the asym- 
metric approximation (21). 

In order to demonstrate the connection between (21) and (77), the latter result must be specialized to an axisym- 
metric mean flow as done by Goldstein [3]. Introducing the mean-flow profiles M(r) and a(r) into (78) shows that 
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Figure 1 0. As figure 8 but for St = 0.5. 


rs (< ^ ) is constant along each ray. If (70) and (7 1) are then used to eliminate 0oo in favor of 0 in the remaining equations, 
they become, to the required order of accuracy, 

r = ± [Q(r \u t ) + cos 2 6] , r 2 <p = r/*, 

S ± = v*<p + <? 2 (r|t/*) - sin 2 0, 

which must be solved subject to 

r = r s , tp = <p s , S ± = 0, 
at r = 0, where = r s q s sin(A — <p s ) = rs M , 

rQ{r |i/.) = s/r 2 q 2 - v*, (33) 

a dot denotes differentiation with respect to the ray parameter r and use has been made of (75) in arriving at the 
definition of rv The order 1 /R term in the f equation must be retained because it leads to an order-one contribution 
in the far-field behavior of S ± . The branch cuts of Q are specified below. For the present, attention will be restricted 
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Figure 11. As figure 8 but for St = 0.25. 


to the case where v+ is purely real which, it turns out, implies Q 2 ^ 0 along the rays. A necessary (but not sufficient) 
condition for to be purely real is that the source be located such that q 2 s > 0. 

The choice of sign in the r equation of (32) is a function of the initial condition, 

r = q s cos(A — <p s ) at r = 0, 

as well as the number of zeroes of Q encountered along a ray trajectory r(r). It is assumed here (in agreement with 
the turning-point assumptions made in Appendices A and B) that rQ is a monotonically increasing function of r with 
at most one simple zero. rQ will then be non-zero and r will increase monotonically with r for rays initially directed 
away from the jet centerline, i.e. -n/2 < A — <p s < tt/2. These solutions, for which the positive sign in (32) applies, 
will be referred to as direct rays. For rays initially directed toward the jet centerline, i.e. n/2 < A — <p s < Sn/2, but 
which eventually reach the far field, r initially decreases with increasing r until r = r$ where Q vanishes. Once this 
occurs, r begins to increase with r and the sign in (32) must change from negative to positive. These solutions will be 
referred to as indirect rays. 
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In view of the preceding discussion, the solutions to (32) can be written as 




v+dr 


r 2 <5(r|i/*) ! 


where 


S ± = v*A(p + £(r|i/,) =F C( r s|^*) ~ -Rsin 2 9 - R x + R , 
C(r|j/*)= f Q(r\u a )dr, Q 2 (r 5 |t/*) = 0, 

J rx 


(34) 

(35) 

(36) 


and the upper (lower) signs in (34) and (35) apply to the direct (indirect) rays. Notice that r$ cancels out of both the 
A<p and S ± solutions when the upper set of signs are chosen. Substituting these results into (77) yields 




±r s Q{r s \u i ,)^- 


"I 


W * ) — Rsin 2 6] 


(37) 


as fc 0 , R — ► oo, where d^* = ±r 5 Q(r s |i/*)dA follows from the definition of v+. 

The connection between (37) and (21) is established by showing that the former result is nothing more than the 
leading-order, large - ko approximation of the latter. This is done by using the Poisson sum formula [14] to rewrite (21) 
as 


Ku. 


+°° /»-j-oo 
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~ J — oo 
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r s Q{r s \v) 


Ai[r ? (r s |^)]e iA:o[Q ’" 1 ' +c(r|l/)_ ' Rsm2 e] du, 


(38) 


as ko, R oo, where a m = A <p 4- 27 rm, 


V(r\v) = ~ [l^oCMAI 3 » 


and Q{r\v) and Q{r\u) are given by (33) and (36) with replaced by v. 

Since the argument of the square root in (33) can now be negative (or even complex when the v integration of 
(38) is performed in the complex plane), the branch cuts of Q must be made explicit. In the present context, it is clear 
that Q is simply a generalization of the function Q n defined by (16) where the square root was chosen so that it has a 
positive imaginary part for negative (real) arguments. The appropriate generalization of this choice to arbitrary values 
of v is 

rQ{r\u) = y/\ r 2 q 2 - v 2 \ e i iI ar s( r 9-*')+ ax 
with — |7r < arg(r<? ± v) < 

The branch cuts of Q(r s \v) in the complex u plane are shown in figure 12. It is worth noting here that 

V~v( r sW) 


(39) 


lim ^ , . . 

v'-rlql r s Q(r s \u) 


= [ko?(,r a )]* , 


where (P(r) = l/r(r 2 q 2 )', which shows that the integrand in (38) remains bounded at the branch points of Q{r s \v). 

The integral in (38) is evaluated asymptotically in the limit ko — > ► oo using the method of steepest descents [11] 
which requires making the ko dependence of the integrand explicit. The Airy function Ai can be replaced with its 
large argument behavior along the entire real v axis except in the small order k^ 2/3 neighborhoods of the branch 
points ±r s q s . Thus 

' Ry - my +Tiy for q 2 > 0 , 

ny for q 2 < 0, 
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Figure 12. Branch cuts of Q(r s \v) in the complex v plane, (a) q 2 s > 0, C {+) lies between ±r s q 3 where it overlaps C<~); ( b ) Q 2 s< o. 


and 7£^ x) is given by the right-hand side of (38) but with the integration done over the contour C (x) . The contours 
C ( ~\ C (+) and C (x) are shown in figure 12. 

For the present, it will be supposed that the dominant behavior of (38) is not determined by 7£^ x) . This exceptional 
case is considered in the following section. Applying the method of steepest descents to (41) yields 




rsQ(rsK) 


dip 




e i(/3+ 1 n) G ik 0 [a mit +C(r )=fC (r a W*)~ Hsin 2 i 


(42) 


as /co — > oc, where 

-|tt < 0 = - arg 

m* and v+ are determined by the saddle-point condition 


otm, = Ay) + 2 Tim* = V'K) - ( / =F J ) ^2^ 


dr 


2 Q(r\ v *y 


(43) 


and the ★ subscript is reused in order to emphasize the connection with the ray-theory solution. 

It is immediately evident that the ray solution (34) and the saddle-point condition (43) are merely different versions 
of the same relation. The latter result determines v+ as a, possibly multi-valued, function of A <p. Multiple solutions 
for at a fixed A<p are accounted for by the summation in (42) and indicate different rays reaching the same far-field 
observation point. The ray solution (34) determines A <p as a single-valued function of v+ so no special treatment is 
needed for rays that reach the same far-field position. The 2n m* factor in (43) allows for the possibility that 'ip falls 
outside the range [ — 7r, -Hr]. 
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Figure 13. Contours in the complex r plane (a, c ) and corresponding regions in the complex v plane (b, c/). Shaded areas indicate 
location of branch cuts of Q(r|i/) along indicated r contours, (a, b), pertain to £(r|i/*) + C( r sl 1 '*) when Qs > 0; (c, d), pertain to 

CM*'*) - CO'sl*'*) when q % < 0. 


Equations (34) and (43) also imply that dip/dv* = dip/dv+. If v± is purely real as assumed when deriving (37), 
ip is also purely real and the amplitude factor in (42) can be rewritten as 
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It then follows that 7v^7 : corresponds to the direct-ray solution given by the upper signs in (37) and the indirect-ray 
solution will correspond to — i provided 


~r s Q(r 
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The above condition removes the ambiguity in the phase of (37) when the lower signs are taken. The — 7t/2 phase shift 
is a consequence of the indirect ray having passed through the caustic at r = r$ before reaching the far field. The shift 
is left undetermined in the ray-theory solution described in Appendix C and, in general, must be obtained through a 
local analysis near the caustic [15]. 

The above results clearly establish the connection between the asymmetric, high-frequency solution (21) and 
the ray-theory solution (37) when v* is purely real. Since the location in the complex u plane of the saddle point 
determined by (43) is not restricted (other than as required by the method of steepest descents), (40) and (42) show 
how (37) can be generalized to complex rays. The technique of applying the method of steepest descents to a classical 
high-frequency solution in order to guide the generalization of a ray-theory result to complex rays is well known 
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[16, 17]. Less well known are the methods for developing a complex ray theory when no classical high-frequency 
solution is available as, for example, in the case of the uni -directional, transversely sheared, mean flow considered in 
Appendix C. A survey of available approaches for doing just this is given by Chapman et al [18]. 


6 Evaluation of ray-theory solution 

In arriving at the lar ge-/c 0 approximation of TZ^ given by (42), it was implicitly assumed that the integrand in 
(41) can be analytically continued into the complex v plane so as to allow integration along a contour that is (at least 
locally) coincident with the steepest descent paths intersecting at When i/* is complex, care must be taken during 
the evaluation of (42) because r$ is then also complex and the integration of Q in the definition of £ must be done 
along a contour in the complex r plane. For each point along that contour, Q has branch cuts in the v plane determined 
by (39). The r-integration contour must therefore be chosen such that these branch cuts leave a region of analyticity 
in the complex v plane that contains both v+ and the real axis. Figure 13 shows example contours in the complex r 
plane used for the evaluation of C( r W*) + C( r s W*) and C( r W*) — C( r s\v*) as well as the corresponding regions of 
analyticity in the complex v plane. 

Figure 14 contains a plot of the saddle point v+ determined from (43) as a function of polar angle 9 for A <p = 60, 
r s = 0.75 and the mean flow given by (25) and (27) with Mj — 0.9. The direct (indirect) ray solutions, which 
correspond to the upper (lower) signs in (43), are indicated by a A (□). The real and imaginary parts of are denoted 
by the open and closed symbols, respectively. 

The figure shows that the complex solutions for v+ are confined to the range 6 < 39.5. These solutions describe the 
so-called zone of silence where a significant reduction in the sound radiated to the far field results from an exponential 
decay in 7 Z u . Notice that the boundary of the zone of silence does not coincide with the value 9 % 30.67 where 
q 2 s = 0 (cf. figure 2b). Just outside the zone of silence is a range of 9 where multiple (real) solutions for v+ are found. 
These multiple solutions occur when different rays reach the same far-field observation point and give rise to ray- 
interference effects which can be either constructive or destructive depending on the relative phases of the solutions. 
The ray-interference region extends to 9 « 58 beyond which q 2 has no turning points (cf. figure 2b) and only direct-ray 
solutions are found. 

Also shown in figure 14 is the location of the branch point r s q 3 along the real v axis. The real part of z'* intersects 
this curve at 6 « 37.5 with the result that the solution type changes from a complex direct ray ( 9 < 37.5) to a complex 
indirect ray (9 > 37.5). The change in solution type is a consequence of the phase shift in Q{r s \v+) that occurs as 
v* crosses the branch cut issuing from r s q s (cf. figure 12a). The asymptotic approximation (42) remains valid for z/* 
arbitrarily close to the branch cut (but outside the order k $ 2/3 neighborhood of r s q s ) because the contribution to the 
large- /co behavior of (41 ) obtained by deforming the ^-integration contour around the branch point is of higher order. 

The purely real solutions for v+ are also affected by an encounter with the v s q s curve with the result that the 
indirect-ray solution (6 < 42) changes to a direct-ray solution (9 > 42). The change in solution type occurs as rs and r s 
become coincident and corresponds to the change that occurs in (34) as A — (p s passes through n/2. The approximation 
of 7Z^ given by (42) becomes invalid near this point because moves into the order 2/3 neighborhood of a Q(r s |z^* ) 
branch point. The dominant behavior of (38) is then determined by 7£<j° rather than 7^Jj _) or 7Z£K The appropriate 
asymptotic behavior of TZ^ is found by applying the method of steepest descents to (38), where the integration is done 
over the contour C (x) , with the result that 
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Figure 14 reveals another change in solution type which also leads to a local breakdown in the large- ko approx- 
imation given by (42). This occurs at the zone of silence boundary 0 ~ 39.5 where the complex indirect-ray saddle 
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Figure 14. Saddle points for an isothermal jet with r s — 0.75 and £<p = 60. Open symbols, Re*/*; closed symbols, Im^*; A, 
direct-ray solution; □, indirect-ray solution; Solid line, R er s q 3 . 


point changes into a pair of real indirect-ray solutions. At the point of bifurcation, z/* becomes a saddle point of 
higher-order for which 



(47) 


with the consequence that the right-hand side of (42) becomes unbounded. Equation (47) implies a zero in the Jacobian 
determinant J introduced in Appendix C and therefore the appearance of a caustic [3], i.e. an envelope of real ray 
trajectories. The breakdown is restricted to indirect-ray solutions because they are the only ones that encounter a 
caustic before reaching the far field. 

An expression for dxjj/dv* can be derived by using (39) and the defining equation for rs to show that 
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Figure 15. |Gu, /&, | x 10 4 for an isothermal jet with r s = 0.75, Atp = 60 and St = 2. Solid line, asymmetric approximation; dotted 
line, ray-theory solution; dot-dashed line, near-branch-point solution; dot-dot-dot-dashed line, near-caustic solution. 




where (P(r) = 1 /r(r 2 <j 2 )'. Using these results when differentiating (43) with respect to leads to 


/ f°° f r ‘\ 1 + 2 vjr&jr) 2i/ 2 3 > (r 8 ) 

& v * \Jr s ^ Jr 5 ) r 2 <5(r|i/.) r ± r s Q(r s \u.)' 


(48) 


where the singularity at r = rs is integrable. 

A large-Aro approximation of can be constructed when the saddle point approaches a zero of dip/dis* by 
using the procedure outlined in appendix G of reference [19]. The integrand is expanded about the midpoint v rather 
than the saddle point when applying the method of steepest descents where v is determined by 
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Notice that, in view of (48), v is independent of A <p. The resulting asymptotic approximation is 
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Figure 16. Various quantities evaluated at the zone of silence boundary for an isothermal jet with r s - 0.75. ( b ) □, midpoint v\ solid 
line, Ker s qs- 


rj = -fcofy'O') - 


'k 0 

d 2, tp 

2 

dp 2 


and ${&) is given by (43) with zv* replaced by z/. 

Figure 1 5 is a plot of | as a function of polar angle 9 for Strouhal number St = 2 and = 60. The solid 

line corresponds to the asymmetric high-frequency approximadon (21). The dotted line corresponds to the ray-theory 
solution given by (40) and (43). Results based on the near-branch-point solution (44) and the near-caustic solution 
(49) are indicated by the dot-dashed and dot-dot-dot-dashed lines, respectively. The curves were computed using the 
same mean flow and point source location as figure 14. 

Comparing the dotted and solid curves shows that the ray solution is in good agreement with the asymmetric 
approximation over most of the 9 range including 9 < 30 which is well inside the zone of silence where the rays are 
complex. The discrepancy near 9 = 58 is most likely due to the failure of the ray-theory solution to correctly describe 
the disappearance of the indirect-ray contribution (cf. figure 14). The ray-theory result could probably be improved 
by constructing a local solution that accounts for the higher-order zero in Q that emerges for the indirect-ray solution 
near 6 — 58. A similar explanation likely applies for the discrepancy near 9 — 150 since an additional indirect-ray 
solution (with m* = — 1) appears in the range 9 > 150. 

Figure 15 also shows that the ray solution breaks down at the zone of silence boundary 9 ~ 39.5 as expected. 
It is interesting to note that (42) remains bounded (although not in particularly good agreement with the asymmetric 
approximation) at 9 « 42 where v+ equals the branch point value r s q 3 because, as can be shown from (48), 


lim 


r s Q{r s =±2r^?(r,), 


where again the upper (lower) sign corresponds to the direct (indirect) ray solution. 

Comparing the dot-dashed, dot-dot-dot-dashed and solid curves shows that the local solutions (44) and (49) bring 
the ray-theory result into closer agreement with the asymmetric approximation and hence the exact solution in their 
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respective regions of applicability. In the range 35 <9 < 50, the near-caustic solution (dot-dot-dot-dashed line) yields 
better agreement with the asymmetric approximation than does the near-branch-point solution. One might then expect 
that a composite solution formed from (40) and (43) together with the near-caustic solution (49) should produce good 
agreement with the asymmetric approximation over the ranges of 9 and A <p of interest. However, as A(p decreases, 
the near-branch-point solution becomes the better approximation in the vicinity of the zone of silence boundary and a 
different approach to constructing a composite ray-theory solution will be taken here. 

Figure 16(a) contains a plot of 9 as a function of A where 9 is defined as the value of the polar angle 9 at which 
the indirect-ray saddle point satisfies both (43) and (47), i.e. the value at which = v. The 0 curve therefore marks 
the zone of silence boundary in 9- Aip space. As A ip — > 0, 9 approaches 30.67 which is the value of 9 at which q 2 s 
vanishes. As A <p increases, 9 also increases indicating that the zone of silence grows as the azimuthal angle between 
the far-field observation point and the source point increases. 

Figure 16(b) shows plots of the indirect-ray saddle point i/* (which is equivalent to v here) and the location of the 
branch point r s q s along the real v axis as functions of A at the zone of silence boundary. It is clear from the figure 
that v eventually moves into the order /cq 2/3 neighborhood of a Q(r 5 |i/*) branch point as A <p becomes sufficiently 
small. When this occurs, the near-caustic approximation (49) must be reworked in order to account for the presence of 
both a branch point and a higher-order saddle point. Rather than deriving yet another local approximation for 7Z the 
approach taken here is to modify the near-branch-point approximation (44) by shifting the solution to (45) as follows 


v = i>- v \ e=9 + sign(/>) r s q s \ e= g. 

Thus v coincides with the Q(r s |i/*) branch point when 9 = 9. If the integrand in (38) is expanded about 0 rather than 
i> when the method of steepest descents is applied the following expression is obtained 


7Z ix) 

/v u ; 


[right-hand side of (44)] x exp 


. ko (Aip — ^) 2 
2 dip/dv 


(50) 


as ko — ► oo, where 'ip is given by (46) and all occurrences of v must be replaced by v in the above result. 

The dot-dashed line in part (b) of figures 3-6 corresponds to a composite ray solution for \G u; /G u , \ at Aip — 30 
formed from outer and inner expansions given by (40) and (50), respectively. A multiplicative composite form was 
used when r)(r s \0) > 0 and an additive form otherwise [20]. The figures show that the composite ray solution does an 
adequate job of extending (40) through both the branch point and caustic at all St. Apart from the discrepancies near 
9 — 58 and 150 discussed above, the composite solution is in fairly good agreement with the asymmetric approxima- 
tion (and hence the exact result) down to St = 0.5. The disagreement near 9 — 58 and 150 tends to spread over a 
wider 9 range as St decreases which supports the conjecture that a local large- fc 0 solution could improve the ray-theory 
result in these regions. It is interesting to note, however, that the composite ray solution gives a better prediction of 
the exact result near 9 = 58 as St decreases which may indicate that the asymmetric approximation over emphasizes 
the ray interference effects when St is small. At St = 0.25, the composite ray solution continues to accurately predict 
the asymmetric approximation inside the zone of silence, but, outside the zone of silence, it predicts a more gradual 
decline in \G u; /G U ;\ with 0 which puts the ray-theory result in better agreement with the exact solution. 

A determination of the level of agreement between the exact and high-frequency asymptotic Lilley’s equation 
Green’s function for more general mean flows would require using the methods of reference [18] to extend the analysis 
of Appendix C to complex rays and then supplementing those results with local solutions of the sort described in 
reference [15] near any caustics or branch points. Such a program will not be undertaken here. Instead, it may be 
inferred from the success of the composite ray solution at predicting the asymmetric approximation that the level of 
agreement between the exact and high-frequency asymptotic solutions for more general mean flows would be similar 
to that shown in figures 3-1 1 between the exact solution and the asymmetric high-frequency approximation. 


7 Conclusions 

It has been shown that, for parallel round jets, the asymmetric high-frequency approximation, which applies to 
sources away from the jet axis, provides the best overall prediction of the exact Lilley’s equation Green's function and 
remains accurate for Strouhal numbers as small as 1/2. The quasi-symmetric high-frequency approximation, which 
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is arrived at by making a near-axis source assumption, was found to be most successful when applied to the ring- 
source directivity and gives a good approximation of the exact ring-source result at all Strouhal numbers considered 
provided the polar angle from the downstream axis is sufficiently large. In this range of angles, an equivalence of the 
quasi-symmetric and asymmetric ring-source approximations was demonstrated. 

The ray-theory solution was shown to be closely connected to the asymmetric high-frequency approximation and 
this close association was used to guide a generalization of the ray-theory result to complex rays. When combined 
with appropriate local solutions near the caustic and branch point, the ray-theory solution was found to be in good 
agreement with the asymmetric approximation and hence the exact result for the Lilley’s equation Green’s function. 
This finding was used to infer the potential for success of the high-frequency asymptotic Green’s function for more 
general (i.e. non-axisymmetric and/or non-parallel) mean flows. 

The numerical results presented here apply to stationary sources embedded in a subsonic parallel round jet and are 
restricted to flow situations where Lilley's equation contains at most one simple turning point. The generalization to 
multiple and/or higher-order turning points would require modifying the WKB analysis given in Appendix A but this 
is straight-forward and presents no great difficulty other than algebraic. The only impediment to applying the results 
of the present analysis to supersonic flows is the possibility of encountering a zero in the denominator of (24) when 
the source is located such that M s > 1. This singularity can however be ‘removed* by the techniques developed by 
Ffowcs Williams [21]. Extension to sources convecting in the mean-flow direction can be made by simply introducing 
an appropriate Galilean transform into (4). 

It was noted that development of a uniformly valid high-frequency asymptotic Green's function for more general 
mean flows would involve extending the ray-theory analysis of Appendix C to complex rays and then supplementing 
those results with appropriate local solutions near any caustics or branch points. Although mathematically possible, 
the result of such an approach may not lend itself well to implementation in a jet-noise prediction scheme because the 
Green's function would then be given as a function of the initial ray direction rather than the orientation of the far-field 
observation point. One possible resolution to this difficulty may be the direct numerical integration of the Eikonal 
equation (61) and the amplitude equation (62). 
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A Asymmetric, high-frequency approximation 

In this appendix, high-frequency asymptotic solutions for v\ and V 2 are constructed using the azimuthal wavenum- 
ber scaling n = 0(kf). The solution forms are strongly dependent on the turning points of (13) which are determined 
by the zeroes of Q\. It follows from (16) that 


{ —(n/kor) 2 as r — » 0. 
sin 2 6 as r — > oo, 


and therefore that (13) always has at least one turning point with the present scaling of n. 

For isothermal, subsonic jets with monotonically decreasing Mach number profiles, it can be shown [6] that Q\ 
has at most one simple zero for 0 < 6 ^ 90. The situation becomes complicated for more general jet profiles and for 
polar angles in the range 90 < & < 180 due to the possibility of multiple and/or higher-order turning points which 
then depend on the detailed shape of the mean-flow profiles. In order to keep the analysis as simple as possible, it will 
be assumed here (as well as in Appendix B) that (13) has at most one simple turning point. The present analysis can, 
if necessary, be extended to more complicated situations by making some straight-forward modifications to the results 
given here. 

Let r$ denote the single n-dependent turning point of ( 1 3) then the general solution 


v ~ < 


_ 2 

r 6 - r » k 0 3 , 


|<2 n |-l(A-e- i/coC " +^-e ifco ^), 
lAi(7?„) + BBi{rj n ), \r - r*| = Oik^ i ), 


(51) 


|Q n |-4(A + e ifco ^ +B + e-' ,ko <"), 

as /cq — ^ oo, follows from WKB theory [11], where 


V n = x(rs-r), x= Uo 


r - rs » k { 


0 > 


r=7\5, 


C n(r) = [ Qn(r)dr, Q 2 n (r s ) = 0, 

Jr s 


and Ai and Bi denote Airy functions in the notation of reference [12]. The constants A ± , B ± , A and B are related by 

h 


2 A~ = 2 A = A+ji” + B + e~''**, 


(52) 
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B~ = 


(53) 


- 7 -) 2 B = A + e~ i i n + B + e‘i\ 

7T/C 0 J 

which ensure matching between the limiting forms in (51). 

Applying the outgoing- wave condition (14) to (51) shows that B+ = 0. The corresponding solution can then be 
written as 

vi (r) ~ Q n ^ (r)A^ 

as ko — > oo, where attention is restricted to the r - r$ » 2/3 behavior since that is all that is required in (20). 

The bounded condition (15) requires that B~ =0. In this case, only the value of the corresponding solution at 
r = r s appears in (20). However, the location of the turning point r$ relative to r s varies with both 0 and n and it is 
therefore convenient to express the V 2 solution in the uniformly valid composite form 


v 2 (r s ) ~ A 2 


Qn( r s ) 


Ai[? 7 n (r 5 )]. 


as ko —> oo, where 


Vn (r) = -[lkoUr)] 2 \ 


with the cube root defined such that 7] n ^ 0 for r ^r 6 . 

Since the Wronskian V is independent of r, the r — r$ k 0 2/3 behavior of i>i and v 2 can be used to show that 

V = Viv' 2 — v[v 2 ~ — YlkoA'l A 2 e 1 ^, 

as ko —i > oo. Combining this result with the expressions for V] (r) and v 2 (r s ) given above leads to 


v\ (r)v 2 (r s ) 

V 


i7T y/-?7n(rs) 
*0 Qn{'T's')Qn(l'') 


Ai [r)n(r s )]e 


i^oCn(r) 


(54) 


as fco — > oo. 


B Quasi-symmetric, high-frequency approximation 

In this appendix, the high-frequency solutions to (13) are constructed using the azimuthal wavenumber scaling 
n = 0(1). Equations (16) and (17) show that k^Q^ and § then become of equal order as r becomes sufficiently small. 
The disordering of (13) is dealt with, as in reference [2], by introducing an inner region where 


r = k 0 r = 0(1). 


The corresponding equation for v is given to the required order of accuracy by 

f2 ^ + [ f V( °) - n 2 + i] v = 0, 

which has the general solution 

v~Vf[AJ n (0 + BH^(i)}, (55) 

as ko — » oo, where £ — q(0)r and J„ and H^ 1} denote the Bessel and Hankel functions of the first kind, respectively. 
When r is order one, the solution to (13) is determined, to the required order of accuracy, by 

v" + klq 2 v = 0, (56) 

which can, of course, be solved using WKB theory [11]. The particular form of the solution depends on the number 
and nature of the turning points determined by the zeroes of q 2 . As in Appendix A, it will be assumed here, for 
simplicity, that the v equation has at most one simple turning point. 
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First, suppose that (56) has a turning point at r = r a > 1/ko where the notation r G is used to distinguish the 
zeroes of q 2 from the n-dependent zeroes of Q\. The solution to (56) is then given by (51), (52) and (53) with n set 
equal to zero. Matching that result with (55) as f — > O(k 0 ) (noting that f is purely imaginary) requires 

\PhxA~ = e' i(fco ^“i n7r) i, y/2 nB~ = -i2e i(/c °^“i n7r) £ ? 


where 


$(r) = / q(r)dr, 
Jo 


and a subscript a is used to indicate evaluation at r — r a . 

Applying the outgoing-wave condition (14) shows that B+ — 0 from which follows 


«i(r) ~ A+ ( 2 H<, 1 >(* 0 Oe~ ,( * o{ '~ 4T,,r) , 


(57) 


as ko — > ► oc, where attention is restricted to the r — r a k ^ 2/3 behavior since that is all that is required in (20) and 

the Hankel function has been introduced (without loss of generality) in order to facilitate the application of Graf's 
addition theorem [12] in section 3. 

The bounded condition (15) requires B = 0 and the corresponding expression for v 2 (r s ) depends on the location 
of r a relative to r s . Reasoning that the quasi -symmetric approximation is only appropriate when r s <C 1, Balsa [2] 
employed the simplifying assumption that r a —r s > 2/3 , i.e. that the source is always closer to the jet axis than the 

turning point. Using this assumption, v 2 (r s ) can be expressed as 


v 2 (r s ) ~ A 2 



i 


(58) 


as ko — » ► 00 . It should be noted that for any given value of r s there will be, in general, a range of 6 for which r s > r a 
and, consequently, for which the above expression is invalid. The implications of this failure are discussed in section 
4 where comparisons with the order-one frequency solution for the reduced Green’s function are given. 

Unlike the situation encountered in Appendix A, it is possible, with the present scaling of n, that the equation 
governing v has no turning points. When this is the case, the solution to (56) is given as 

v ~ q-% (A + e lfeo<0 + B + e" ifco<0 ), 

as > 00 . Matching with (55) as f — * O(ko) (noting that £ is now purely real) requires 

Vi2^A+ = (A + 2B)e~‘i nn , \/i2 nB + = i Ae^ nn . 


Expressions for v\ (r) and t >2 (r s ) can be derived by applying the boundary conditions (14) and (15) as done above. 
It turns out that the final expressions are in exact agreement with those given by (57) and (58) if the convention that 
r a = 0 when q 2 has no zeroes is adopted. 

The r-independent Wronskian V can be evaluated using the r — r a > k^ 2/3 behaviors of v\ and v 2 and is given, 
for the one- and no-turning point solutions, by 

V = vW 2 - v[v 2 ~ -ifco i AtA 2 e- i( ~ ko ^-i nn \ 
as ko — ► oc. Combining this result with the expressions for v\ (r) and v 2 (r 3 ) given above leads to 

V - lir ¥ M ~t(~Y (»> 

v 2 \QsQj 


as ko — ' ► 00 . 
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C Ray-theory approximation 

In this appendix, the high-frequency solution to (4) obtained from ray theory is reviewed. Following the matched 
asymptotic analysis given by Durbin [4], the solution to (4) for the uni -directional, transversely sheared mean flow (1) 
is 

Gu,(x\x s ) ~ ^(x)A(x\x s )t lkoS ^ Xs \ (60) 


as ko — > oc, where the Eikenol S satisfies 


the amplitude function A satisfies 


£ 2 


\ s \ = o, 


seV5. 


.M . 

s + z — <P ) A“ 
a 


= 0 . 


(61) 


(62) 


and here = (1 — Mi*s)/a. 

The first-order partial-differential equation (6 1 ) is reduced to the coupled system of ordinary differential equations, 


.M x 7 • 

x — s + 2 — s=±V L (<P). 5 = s-x, 

a z 


(63) 


along the rays x(r) by the method of characteristics, where r is a parameter that varies continuously along the ray, a 
dot indicates differentiation with respect to r and V_l is the gradient operator in the y-z plane. Equations (63) must 
be solved subject to initial conditions at the source, 


x = x s , x = <r s {cos y, sin \jl cos A, sin y sin A}, 5 = 0, 


(64) 


at r = 0, where the free parameters y and A determine the initial ray direction relative to the Cartesian coordinate 
system {x, y , zj and it follows from (61) and (63) that 


Mj sin 2 y. 


(65) 


The amphtude function A is found by solving (62) subject to matching with a near-source solution. It follows from 
the analysis of reference [41 that 


yl(x|x.s) — 


1 


47ra 5 


a'; sin y 


where 


J = 


d(x,y,z) 


( 66 ) 

(67) 


d(r, n, A) 

is the Jacobian determinant. 

When attention is restricted to the far-field behavior, the solution (60) can be simplified by noting that the rays 
x(r) become straight lines in the absence of a mean flow. Thus, in the far field, the approximation, 


x ~ x s + RoolcosOoo^sinOooCos^oo^sinOooSiiKpoo} 


( 68 ) 


can be introduced, where R ^ is the distance between the far-field observation point and the source point and #00 and 
<j ) 00 are the far-field polar and azimuthal angles measured from axes passing through the source point and aligned with 
the x and y directions, respectively. It is important to note that <j> oo, #oc and R ^ are not equal to <p, 9 and R of sections 
2 and 3 but approach these quantities in the far field, i.e. 


$00 

Roc 


r s 

~ + — esc 0 sin A (p -b * • • , 

R 

(69) 

9 + ^cos0 (1 - cos A ip) H , 

R 

(70) 

R -f r 3 sin 0(1 — cos A ip) H , 

(71) 
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( 72 ) 


as R — ► oo. 

It follows from (61) and (63) that R ^ = 1 which can then be used when inserting (68) into (67) to obtain 

T d 2 ■ /I Woo, <Poo) 

J~R 00 sm 9 ao d ^ x) , 

as Roc oo. Since (63) implies that the quantity i*s is constant along each ray, (64) and (68) show that 

a 2 s (j s cos fi — M s 


1*3 


a? - M 2 


= COS^oo, 


(73) 


and furthermore, in view of (65), Oqo = #oo(aO and 


sin Ooc 


d^co 

d/i 


a s a s sin M- 


An additional consequence of (73) is 

a 2 sin 2 ji = - cos 2 Ooo *. 

which follows directly from the Eikenol equation (61). The Eikenol itself has the far-field behavior 

S ~ Roc 4 - S ± (x\x s ), 


(74) 

(75) 

(76) 


as R 0 


oo, where 


= [ {S'X ± - 

Jo 


sin 2 Oqo 


)dr 


remains bounded as Roo — * oo and x ± = {y, z } denotes the ray vector in the y-z plane. 

Substituting (74) into (72) and the result into (66) yields the far-field approximation for A which when substituted, 
together with (76), into (60) leads to (18) and (19) where now 




/ cfy\ 5 ikoiSj.+Roo-R) 

[dXj 


(77) 


as ko, R — > ► oo, and use has been made of the far-field relations (69)-(71). 

For the purposes of the present investigation, it is convenient to restate (63) and (64) in terms of the cylindrical 
coordinates {x, r, <p} of section 2. The ray equations in the y-z plane are then 


(r) 2 — <p 2 ~ cos 2 Oo o — s (v?)2 , np = s M , 

(■ rs <*>) = ±d($ 2 )/dip, 5 ± = rs M ip + (f) 2 - sin 2 Ooo, 

which must be solved subject to 

r = r s , ip — <p s , = cr s sin/isin(A - <p 8 ), S ± — 0, 



at r = 0, where rs M = dS/ dip. 
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